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Abstract 

In this paper we complete the derivations of finite volume partition functions for QCD using 
random matrix theories by calculating the effective low-energy partition function for three- 
dimensional QCD in the adjoint representation from a random matrix theory with the same 
global symmetries. As expected, this case corresponds to Dyson index (5 = 4, that is, the Dirac 
operator can be written in terms of real quaternions. After discussing the issue of defining Ma- 
jorana fermions in Euclidean space, the actual matrix model calculation turns out to be simple. 
We find that the symmetry breaking pattern is 0(2Nf) — > 0(Nf) x 0(Nf), as expected from the 
correspondence between symmetric (super)spaces and random matrix universality classes found 
by Zirnbauer. We also derive the first Leutwyler-Smilga sum rule. 



1 Introduction 



The recent strong interest in the low-energy limit of the QCD spectrum originates in the Banks- 
Casher relation Q 

relating the density of eigenvalues of the Dirac operator at the origin, p(0), to the order parameter 
for spontaneous breaking of chiral symmetry, the quark condensate S, in the thermodynamic and 
chiral limits (V here is the size of the box). It is expected on the basis of lattice QCD simulations 
that chiral symmetry will be restored above some critical temperature and/or chemical potential. 
By studying the Dirac spectrum in the infrared limit, one hopes to obtain analytical information 
about the chiral phase transition. 

The presence of gauge fields alters the distribution of the Dirac eigenvalues. For free quarks in 4d 
the spacing between eigenvalues is ~ V -1 / 4 , whereas for interacting quarks it is ~ V^ 1 . Moreover, 
it is known that the Dirac eigenvalues derived from QCD in a finite volume are constrained by sum 
rules first discovered by Leutwyler and Smilga p. These sum rules are obtained by expanding 
the partition function in powers of the quark mass m before and after averaging over the gauge 
field configurations, and matching powers of m. 

The same sum rules can also be obtained from a random matrix theory with the same global 
symmetries as the QCD partition function. This was first noticed in Q. In the random matrix 
theory (where the size of the random matrices N is identified with the space-time volume V), 
the average over gauge field configurations is substituted by an average over Gaussian distributed 
random matrices with the same structure of matrix elements as the Euclidean space Dirac operator 
"f^D^. The matrix elements can be real, complex or quaternion real, corresponding to Dyson 
indices f3 = 1, 2, 4 respectively. The corresponding matrix models are called the (chiral) Gaussian 
orthogonal, unitary and symplectic ensembles. 

In fact, the extreme infrared limit of the QCD partition function maps onto the same effective 
partition function (usually called the finite volume partition function) as the random matrix 
theory g |, |, |, % g g. In retrospect this is not too surprising, since the chiral Lagrangian to 
lowest order in the quark momenta is completely determined by the pattern of chiral symmetry 
breaking and Lorentz invariance. In the so-called mesoscopic range <C L <C A^, where L 

is the side of the box and A^* is the Compton wavelength of the Goldstone modes, this partition 
function expresses the quark mass dependence in the static limit and in a finite volume. It is 
expressed as an integral over the Goldstone manifold of the composite variables corresponding 
to the pion fields and it is a function of one scaling variable NY,A4, where A4 is a mass matrix 
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in flavor space. In substituting the complicated average over gauge field configurations with 
an average over random matrices, we have the advantage that we can utilize a well-established 
mathematical framework that had its beginnings with the pioneering works in the field by Wigner, 
Dyson and Mehta. 

The sum rules can be expressed using the so-called microscopic spectral density /05(A), defined 
by magnifying the distribution of eigenvalues in the vicinity of the origin (A = 0) on the scale of 



the average eigenvalue spacing. ps(A) is a highly universal quantity ||, ||, [ji], |ll], |T^, 13|. It is 
completely determined by the global symmetries of the partition function, and does not depend 
on the matrix potential, for example. 

However, the sum rules are not sufficient for determining the spectral density itself. Fortunately, 
the Dirac spectrum can be obtained using the so called partially quenched partition function for 



QCD [14 1, in which one introduces bosonic and fermionic valence quark species in addition to the 
ordinary quarks. When the bosonic and fermionic valence quark masses coincide the partially 
quenched partition function coincides with the original QCD partition function. For a nice review 



of this approach see also [15] 



The cases corresponding to fundamental fermions with gauge group SU(2), fundamental fermions 
with gauge group SU(N C ) (N c > 3), and adjoint fermions with gauge group SU(N C ) (N c > 2) 
(labeled respectively by (3 = 1, 2, 4) in four dimensions, and the corresponding cases labeled by 
= 1, 2 in three dimensions have been analyzed in [||, 0, [J. In this paper we will treat the only 
"missing" case, namely three-dimensional QCD in the adjoint representation. We take the color 
gauge group to be SU(N C ) with N c > 2. This case will correspond to a matrix model labeled by 
(3 = 4, as we will find that the Dirac operator is quaternion real, that is, that its matrix elements 
can be written in the form Qfc/ = a?;l + ia^i ■ 3 where the a!ti are real numbers, 1 is the 2x2 unit 
matrix and 3 the triplet of Pauli matrices. Like for all the other cases, we will propose a Gaussian 
random matrix theory corresponding to the low-energy partition function by substituting the 
integral over gauge fields with an integral over a random, quaternion real matrix. By expressing 
the fermion determinant as an integral over Grassmann variables (utilizing the supersymmetric 
formalism developed in fl6|), we will then perform the Gaussian integration corresponding in the 
field theory to the average over gauge fields. After manipulating the partition function further 
and eventually performing the Grassmann integration, we will obtain the finite volume partition 
function. From this partition function the pattern of spontaneous flavor symmetry breaking 
(assuming that such breaking takes place) will emerge. We will assume from the outset that 
there is a nonzero condensate S. The discrete parity symmetry discussed below and in ||. [l7| 
remains unbroken. (For a more detailed discussion of this issue see these references.) 
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2 The parity-invariant Dirac operator in three dimensions 



To begin, let us consider the Minkowski space Lagrangian for QCD in the adjoint representation 

C = --trF 2 + Y,i>f{iy>-rn f )xl> } (2.2) 

/=! 

where F is the gauge field tensor, f) = j^D^, is the covariant derivative for the adjoint 
representation given explicitly in (|4.12j ), and is the quark mass corresponding to flavor /. For 
reasons that will become evident, we call the total number of flavors 2Nf. ipf are quark spinors 
in the adjoint representation and / is the flavor index (the indices corresponding to color and 
spin are suppressed). 

The lowest-dimensional representation of 7^ is given by the Pauli matrices 7 = 03, 7 1 = io\, 
7 2 = ia 2 - In this 2d representation, there is no chiral symmetry, since there is no 2 x 2 matrix 
that anticommutes with the cr^. 

For zero masses m/, the above Lagrangian is invariant under parity P, defined in three dimensions 
by 



tp(t,xi,x 2 ) 
A (t,xi,x 2 ) 
A 1 (t,x 1 ,x 2 ) 
A 2 (t,xi,x 2 ) 



-> j 1 ip(t,-xi,x 2 ) 

-> A (t,-X 1 ,x 2 ) 

-> -A 1 (t,-xi,x 2 ) 

-> A 2 (t,-x 1 ,x 2 ) 



The mass term breaks this P invariance. However, by choosing half of the masses equal to +m 
and half equal to —m, we can achieve a (P, Z2)-invariant Lagrangian || 18]. In terms of 2— 
spinors, this choice corresponds to Nf 2-spinors <j)f with mass +m, and Nf 2-spinors Xf with 
mass —m. Under P the mass terms for the 2-spinors change sign, so that if the two sets of 
two-spinors transform into each other in a Z 2 transformation <j)f Xfi f = 1, 2, ...,Nt, the total 
Lagrangian is invariant under the combined transformations P and Z 2 . We can use this choice 
of mass term to write down a (P, Z2)-invariant Lagrangian in the adjoint representation: 
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1 2N f _ N f _ 2N f 

C = - jtrF 2 + ^2 ^pP4>f ~ m ^f^f + mtpfTpf (2.4) 
/=1 /=1 /=A f /+l 

As we will see in section |], the given representation of - in this case the adjoint representation 
- uniquely defines (in Euclidean space) the anti-unitary operator Q that commutes with the Dirac 
operator il/). We will see that for the adjoint representation, the condition Qip = ip leads in the 
fermionic partition function to an integral over only half of the number of fermionic degrees of 
freedom with respect to the fundamental representation. The condition Qip = ip is therefore 
called the Majorana condition. 

Also for (3 = 1 (fundamental fermions with SU{2) color symmetry) we have an anti-unitary 
symmetry [if),Q] = (with the difference that Q 2 = +1, whereas in the adjoint case we have 
Q 2 = —1; see below), which leads to the same kind of relation Qip = ip. In the case of the 
fundamental SU (2) representation this relation defines a basis in which the Dirac operator has 
real matrix elements. In the adjoint case, because Q 2 = —1, we will find that we can write the 
Dirac operator in terms of real quaternions (see section |5|). In contrast to the adjoint case, in 
the fundamental SU(2) case we integrate in the partition function over both ip and ip, which are 
independent degrees of freedom. This was shown explicitly in section 3 of ref. ||]. For the gauge 
groups SU(N C ) (N c > 3) and fundamental fermions, we have no anti-unitary symmetry. 



To discuss the theory (2.4) in Euclidean space, we first discuss the general issue of defining 



Majorana fermions in Euclidean space. 



3 Defining Majorana fermions in Euclidean space 



It is not in general straightforward to define Majorana fermions in Euclidean space. In 4d, this is 
due to the fact that the would-be Majorana fermions do not transform like Dirac spinors under 
Euclidean Lorentz-transformations. This in turn is due to the absence, in Euclidean space, of an 
equation relating the right- and lefthanded Lorentz-transformations. For Minkowski space there 
is such a relation: 



a 2 A R a 2 = A* L (3.5) 

This relation guarantees, in 4d Minkowski space, that the upper and lower components of the 
usual Majorana fermions, £l and — o"2^2 transform like lefthanded and righthanded components, 
respectively. In || it was shown that one can nevertheless construct a fermionic partition function 
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for Majorana fermions in 4d (or 2d) Euclidean space. To begin with, we will here briefly recall 
how this was done. 

Given the Euclidean Dirac operator for adjoint fermions, ilj) (defined as usual by lp = 7u-D« 
where 7 M are Euclidean gamma matrices satisfying {7^,7^} = 25^ and D^ ab = d^S^ + fabcA'p, 
an anti-unitary operator Q defined by 



[ip,Q]=0 (3.6) 
was identified. Since in 4d it turns out that Q 2 = —1, the Majorana condition 



= i/> (3.7) 
is contradictory (since it implies that tp = —ip and thus ip = 0), unless we define 



ip** = -V> (3.8) 

This condition is called conjugation of the second kind. It is common in the literature on 
Grassmann variables and in calculations involving supersymmetric random matrix theories. If Q 
commutes with Euclidean Lorentz-transformations, a partition function for Euclidean Majorana 
fermions can now be defined if we can find an operator O such that 



Z = y/det(ip) = jD^e-^ Toi ^ (3.9) 

where ip are Majorana fermions satisfying Qip = ip, det O = 1 and Oiip is an antisymmetric 
operator. This last condition guarantees that the square root is well-defined. Lorentz-invariance 
of this partition function is guaranteed because in addition, we demand that O by construction 
satisfies 



^ T Q = ^ = (QV) f (3.10) 

(we recall that in Euclidean space the Lorentz-invariant quantity is tp^ip). In ( |3.9| ) we integrate 
over only half of the degrees of freedom with respect to the partition function for the usual Dirac 
fermions. 
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Conversely, the Majorana condition can be identified by the condition (|3.10| ), once we have written 
down the partition function, by demanding Lorentz-invariance of the latter. 



4 The fermion determinant 



Following these hints, we can now construct a partition function for Euclidean Majorana fermions 
in 3d. As it turns out, the problems with defining Euclidean Majorana fermions that we encounter 



in four dimensions, are not present in 3d because the Majorana fermions defined by (3/7) transform 
correctly under 3d Euclidean Lorentz-transformations. This is because in 3d, there is just one kind 
of 2-spinor and just one two-dimensional representation of the Lorentz group, which is equivalent 
to the rotation group. 

We start by defining Hermitian gamma matrices for our 3d Euclidean space. They can simply be 
taken to be the Pauli matrices 



70 = 03, 7i = o"l 5 72=^2 (4.11) 
satisfying {7^,7^} = 25^ v . The Dirac operator in the adjoint representation is 

l^D^ = J^Sab + iA c ^T c ) ab ) = 7m (0 a A 6 + fabcA^) (4.12) 

where f a b c are real structure constants for the gauge group. is antisymmetric under transpo- 
sition and If) is antihermitean. It is easy to find the antiunitary operator that satisfies ( |3.6| ) in 
this case. It is given by Q = i^iK where 272 = C is the charge conjugation operator satisfying 
C7*C _1 = — 7^ and K denotes complex conjugation. We see that like in the 4d case outlined 
above, Q 2 = — 1 and the Majorana condition 

V> = CKi) (4.13) 

makes sense only if we introduce conjugation of the second kind, ip** = —tp. The Majorana 
condition is then consistent with the explicit form 



H4- i4M) 
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of the spinors. The generators of the Euclidean Lorentz-group corresponding to the representation 
( 8,56] ) of 7 M are given by [19 



S^ = \h„lu] (4.15) 

We find Sq\ = — (T2/2, = <7i/2 and S\2 = — cr 3 /2. These are equivalent to the generators 
of the rotation group and one can easily verify that the corresponding Lorentz-transformations 
commute with the anti-unitary operator Q defining the Majorana condition in our case. 

To write down the fermionic action we now look for an operator O such that det O = 1 and 
ip T = Tp> = (Qip)^ . We immediately arrive at 



O = -ia 2 (4.16) 
In addition, Oil/) is antisymmetric so it is now straightforward to write down a partition function 



for Majorana fermions satisfying ( 4.13 ) in 3d Euclidean space 



det(i#) = y / det(Oi^) = jD^e-^ T ° 1 ^ (4.17) 



5 The Dirac operator for Majorana fermions 



In analogy with 4d, we expect that there is a basis in which ilf) can be written in the form of real 
quaternions |J. We will now verify this. It is easily checked (using C 2 = —1, x*k = ~Xk and 
= that an expansion of ip consistent with the Majorana condition is 



ip = T,(^Xk + C0lxt) (5.18) 



where (fik are arbitrary c- number 2-spinors and Xk are Grassmann variables. Since (CK) 2 = — 1, 
(f>k and C(\)* k are linearly independent. 

We then see that the expression in the fermionic action involving the Dirac operator can be 
written 
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kl ^ ' 

v / Xk Y ( 4>\iP4>i fypctf \ ( xi 



The matrix 



has the form 



(5.19) 



Qkl S ( Sfh $*f C K ) (5.20) 



Qkl = [ % ) (5.21) 



This can be verified using the properties of C and eq. ( |3.6|) from which we derive 



chipc = {ip)* 

tfifl = -(ip)*C (5.22) 



Qki can be written as a real quaternion: 

Qm = a kl l + ia k i ■ a (5.23) 

where a kl , a l kl are real numbers. 

6 Random matrix theory 



We are now ready to define the random matrix theory for massive, parity-invariant three- 
dimensional QCD with fermions in the adjoint representation. In analogy with all the other 



8 



cases in 3d J?], [|] and in 4d || |4|, |5|, |T^] , we substitute the integral over gauge field configurations 
in the partition function for QCD with an integral over a Hermitian random matrix T. We take 
this matrix to have the quaternion structure ( 5.20| ). Using the same notation as in Q we therefore 
set 



a ij ia ^ 

M=0 



(6.24) 



where a^- are real numbers and we have defined ao = —i and o\, 02, 03 are the usual Pauli 
matrices. This means the matrix elements of T are themselves 2x2 matrices. The operator 1J) is 
antihermitian, so we can substitute Ifi in the Euclidean fermion determinant with the matrix iT 
where T is Hermitian: 



Z( mi ,...,m 2Nf ) = f DAe~ s ^ [J ^det^ + m/) 

/=i 

- J DTe~ Ns2t *W J] Jdet(iT + m f ) (6.25) 



This random matrix theory has the same global flavor symmetry as the QCD partition function. 
In addition, T is here taken to be a matrix of N x N real quaternions, so that the anti-unitary 
symmetry of adjoint QCD is reproduced, and as we will see, the flavor symmetry breaking pattern 
will appear from the random matrix partition function. DT is the invariant (Haar) measure. As 
usual I, |, |, g, |, [Uj N is identified with the space-time volume. We will see that E is the value 
of the condensate (the order parameter for spontaneous symmetry breaking) defined by 



1 d 

S = — lim lim — — hxZYmi, .... m2N f ) (6.26) 

m/->0 jV-»oo N omf 1 

We are assuming that this condensate is non-zero. As discussed in Section ||, we choose the masses 
in pairs of opposite sign, so that the corresponding Minkowski space theory is parity-invariant. 
We call the total number of flavors 2iVy, and we choose m/ = +m for / = 1, ...,Nf and mj = —m 
for / = Nf + 1, 2A r j, where m is a real positive number. 



In order to evaluate the partition function ( 6.25| ) we rewrite the square root of the fermion 



determinant as an integral over Grassmann variables (cf. ( |5.19| )): 
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Y[^Jdet(iT + m f ) = j JJlty/exp 

f f 



'E"'/' r/ - <">f )ijo j f 



(6.27) 



where the indices i, j are summed over from 1 to N. Here we are using the formalism developed 
in |l6[] for super symmetric matrix integrals. Our matrix integral is pure fermionic and involves 
no bosonic variables, so we will apply this formalism using only the fermion-fermion block. Our 
integration measure is 



2N f N 

n^/=nn^/V/ 



(6.28) 



f=i i=l 



Recall that Tij are quaternions, so the (jfj are 2-component vectors like in ( |5.19| ): 



4 * 

X) 



(6.29) 



We use conjugation of the second kind x** = ~X f° r the Grassmann variables. 

The first step is to perform the Gaussian integral over the random matrix T in the integral 



Z(m) 



DT f Y[D(p f exp 

f 

Da^ Y[ Dfyexp 
f 



-NY?ix(T*T) -i^4> f *{T- im f )4> f 



- N ^« + E <>'i"W : >f - E V/ 



(6.30) 



where we sum over repeated indices \i and i, j. In the second step of (|6.30|) we have used that 
a^al + <J v cft = l&iiw To perform the Gaussian integral we complete the square in the exponent 



of ( |6.30D by setting 



a ij- a ij 2S 2 7V 



E4> 



(6.31) 
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The symmetry properties of the afj are the same as the symmetry properties of a^, namely 



a-j = a®}, a!?j = — e& (k = 1, 2, 3). These symmetry properties follow from the hermiticity of the 
matrix T and ( |6.24| ). The integral over T is uniformly convergent in the fermionic variables 0^, 
so we can interchange the two integrals and perform the Gaussian integration. We then arrive at 



Z(m) ~ J ^Q-D^yexp 



fg f 



(6.32) 



To evaluate 4> l f $g we use the Fierz' identity 



5>)f al s = 2{8 aS 8 lP -8 a/3 6, 



■ys) 



(6.33) 



and insert the explicit form for the 2-component Grassmann variables 0j 



X) 

i * 



(6.34) 



to arrive at 



(6.35) 



where we have set 



*fg = Xf Xg+ Xg Xf 



(6.36) 



(a sum over i is understood). Disregarding for a moment the mass term, we have rewritten the 
exponent in Z{m) as a square. We can therefore use the Hubbard-Stratonovitch transformation 



exp [-<xF /g F /g ] ~ j da fg exp 



~4^°~ /9<T/9 ~ i(Tf9Ff9 



(6.37) 
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where <Jt g is a real variable, to rewrite the partition function as 



Z(m) ~ / ]^[ Dx/Da exp 
f 



Z 2 N 



a % + a f9 F fg ~ 2 zZ m fX)*x) 
f 



(6.38) 



In ( 6.3Sj ), Da is the Haar measure for the real symmetric matrix a. In order to preserve the flavor 
symmetry of Z(m), a is chosen symmetric like the matrix F. Interchanging again the Grassmann 
integral with the integral over a, and writing out the Grassmann components of Ff g we see that 
the partition function is proportional to 



Z{M)~ jDaUDxfexp — tr (aa T ) + 2 X i > + M) X l 



(6.39) 



where A4 is the 2Nf x 2Nf mass matrix in flavor space, 



/ m 



M 



rn 



-771 



-m J 



(6.40) 



Performing the Grassmann integrals we get 



(6.41) 



with S a symmetric real matrix and DS the Haar measure. In the next section we will evaluate 
this expression using a saddle point analysis to find the low-energy effective partition function 
for three-dimensional QCD with adjoint fermions. 



7 The effective partition function 



We now decompose the real symmetric matrix S in (|6.41|) into "polar" coordinates (2 
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S = OAO T 



(7.42) 



where O is a real orthogonal matrix and A is the real diagonal matrix 



A 



(7.43) 



\2N f ) 



One can always choose the polar coordinates such that 



Ai > A 2 > ... > A2jv f 



(7.44) 



The integration variables in S and those in OAO T will be in one-to-one correspondence if the 
integral over DS is taken to be 



DS= I DODAJ(A) 
oe[o] 



(7.45) 



where [O] denotes the set of left cosets of the group of 2Nf x 2Nf real orthogonal matrices with 
respect to the subgroup consisting of matrices of the form 



/ ±1 



V 



±1 



±1 



(7.46) 



and J (A) is the Jacobian corresponding to ( 7.42| ). This Jacobian was given in and is propor- 
tional to 



J{A)<xJl\Xi-Xj\ 

i<j 



(7.47) 



We are now ready to determine the saddle point of the partition function at zero mass. Setting 
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M = in ( |6.41|) and making the substitution ( |7.42| ), ( |7.45| ) in Z(0) the partition function takes 
the form 



Z(0) ~ f DAexp 



S 2 iV 



(7.48) 



The matrix A is diagonal and the saddle point is given by 



(7.49) 



since the Jacobian drops out in the large N limit. It follows from the derivation of the Banks- 
Casher formula that if the flavor symmetry is broken spontaneously, the condensate for each 
flavor has to have the same sign as the mass Q 0, [|, so we choose, in accordance with (7.44) and 

ED 



/ 1 



Asp — , „ , 



V 



(7.50) 



1 / 



We now expand the determinant in ( |6.41 ) at the saddle point, A = A sp for a small mass matrix 
Ai 7^ 0. To first order in Ai we find 



Z(M) ~ J DOdet N {OA sp O T + M) 

~ J DO&et N {Ok sp O T )e m ^ 1+0 ^° TM ) 
oc J DOe NT,tl( - OJ ° TM) 



(7.51) 



The matrix J was defined in ( ff750|) . This is our expression for the low-energy effective partition 
function. As usual, it is a function of the scaling variable NUM. The matrix J is invariant under 
the subgroup O(Nj) x O(Nf), so in our final expression for Z(M), 
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Z(M)~ DOe NBtv(OJO^M) {7 52) 

J0(2N f )/(0(N f )x0(N f )) 

the integral goes over the coset space 0(2Nf)/(0(Nf) x 0(Nf)). The flavor symmetry breaking 
pattern is thus 0(2Nf) — > 0(Nf) x 0(Nf). This is also what we expect from the fermionic action, 



2N f 



S F = [ d 3 x ^ ip]0{ip + m f )ip f (7.53) 

It is evident that J^f^J® is invariant under 0{2Nj) transformations in flavor space (note 
that the operator Olfi is diagonal in flavor space), while the mass term J2f rrifipjipf 1S invariant 
under 0(Nf) x 0(Nf), which is the unbroken subgroup. The dimension of the coset is 

M= 2iV / (2^-l)_ 2 JV / (iV / -l) =iV/2 (?54) 

8 Sum rules 



We can easily derive the first Leutwyler-Smilga like sum rule for the eigenvalues of the Dirac 
operator. We will use the same method as in || [| . 



The sum rules are obtained by expanding the expression for Z(«M), eq. ( 7.52; ) and comparing the 
coefficients order by order in m 2 to the (normalized) expectation value of the fermion determinant: 



(Note that because of our choice of (P, ./^-invariant Lagrangian, there is an effective chiral 
symmetry that makes the spectrum symmetric around A = 0. This chiral symmetry applies to 
the 4x4 gamma matrices given in eq. (2.4) of ref. II: 



7 1 = 7 J, 7 2 = I T r I (8.56) 
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or equivalently to the corresponding Euclidean gamma matrices. We could just as well write C 
in terms of these.) Here the expectation value is defined as 



/ A ,m } = - Kf f ' ^ ' (8.57) 

jDAe-^(U k jXf)f(\0) 

where A is the gauge field and S[A] the Euclidean Yang-Mills action. Expanding the integrand 
m (gjjp the surviving group integrals at order m 2 have the form 

((X) = f DO tr (O JO T X)ti* (OJO T X) = f DOtr 2 (OJO T X) (8.58) 
JOdG/H JO€G/H 

where G/H is the coset and X = NT,M = iVEmJ. The first order term is killed by the group 
integration. We note that the matrices OJO T are 2Nf x 2Nf symmetric unimodular matrices. 
We now choose real, symmetric and traceless generators t k , k = 1,...,M S for these. M s is the 
number of such generators: 



2Nf(2N f + 1) . 

Ms = f ~ 2 f ' ~ 1 (8.59) 

The generators can be normalized as follows: 

tr(t fc t,) = 2N f 6 kl (8.60) 

This normalization is natural since the generators are 2Nf x 2Nj matrices. We wish to choose 
t\ such that t\ = J. In doing this we can use the relation 



53 tT(At k )tv(Bt k ) = 2N f tv(AB) (8.61) 
fc=i 



to calculate the value of C(^0- Equation (^.61) is valid for any two symmetric unimodular 



matrices A and B. Because of the invariance of the measure in the integral defining (,(X), one 
finds that ((X) is proportional to tr(X'X). Due to eq ( 8.60Q and the choice of t\, we find 
C(t\) = C(^) = ••• = C(*A/)- Therefore we can set 
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C(*l) = TfYLI DOtr\OJO T t k ) (8.62) 
Ms k=i J 

Using dj6g ) and tr(J 2 ) = 2JV) we now immediately see that 

C(X) = ^vol{G/H){NVmf{2N f f (8.63) 
Inserting this into the expansion we get 



' - A + - 2 22V> £ i + -\ = 1 + ^(iVSm) 2 -L(2iV /)2 + ... ( 8 . 64 ) 



Z(0) \ - I 2^ M, 



s 



where the volume of the coset cancels in the ratio. Inserting the value of M s we therefore arrive 
at the sum rule 



\£ mXkf/ " 2(2N f - l)(N f + 1) (8 ' 65) 

where 2Nf is the original number of flavors. Considering also the sum rules found in [0, ||| for 
Dyson index j3 = 2, 1 in 3d, the sum rules for three dimensions can be summarized as follows: 



\5o(^) 2 / 2(2iV7 -1)(^L + i) (8 ' 66) 

where 2Nf is the number of flavors and (3 the Dyson index. The corresponding formula for the 
4d cases was given in [||, ||. In the published version of this paper, as well as in the previous 
electronic version, the wrong sum rule was given due to an error in ( |8.61| ). 



9 Conclusion 



We have obtained the mass dependence of the finite volume partition function for three-dimensional 
QCD with quarks in the adjoint representation and gauge group SU(N C ) (N c > 2). We chose the 
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quark masses such that the corresponding Minkowski space QCD Lagrangian is (P, Z2) -invariant, 
where P denotes the parity transformation defined in three dimensions by eq. (|2.3|) . 

As a starting point we used a random matrix theory with the same global symmetries as the gauge 
theory. This is possible since the gauge theory and the random matrix theory are equivalent in 
the mesoscopic regime. 

The effective partition function describes the static limit of the Goldstone modes resulting from 
the spontaneous breaking of global flavor symmetry. It is determined by the symmetry breaking 
pattern. We assumed that flavor symmetry breaking occurs, and found that the global flavor 
symmetry 0(2Nf) is broken by the vacuum state to 0(Nf) x 0(Nf) symmetry. 



Although the author does not claim to understand the classification of symmetric superspaces, 
it is interesting to note that the result obtained here for the Goldstone manifold, together with 
the other Goldstone manifolds obtained for (5 = 1, 2, and 4 in 3d and 4d, exhaust the scheme of 
Zirnbauer [^l|] with respect to the fermionic symmetric spaces Mp. (I thank Poul Damgaard and 
Jac Verbaarschot for this remark.) We can see this by looking at Table 3 of [21] and matching the 
various random matrix ensembles with the fermionic symmetric spaces given by Zirnbauer. These 
spaces, together with the corresponding bosonic symmetric spaces Mg, make up the Riemannian 
symmetric superspaces corresponding to the various supersymmetric generating functions. In 
simple terms, each symmetric superspace corresponds to a random matrix theory. For us the 
Mp define the respective Goldstone manifolds, since we do not have any bosonic sector. From 
the fermionic symmetric spaces given by Zirnbauer [21] we confirm that our case, the Gaussian 
symplectic ensemble, corresponds to the random matrix theory denoted All in Table 3 of [21], 
while the coset spaces obtained in |?j, || nicely fit the compact symmetric spaces belonging to 
the classes labelled by A, AI, AIII, BDI, and CII. The classes A, AI, and All correspond to the 
non-chiral ensembles while AIII, BDI, and CII correspond to the chiral ones. 
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